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Abstract 

A  quantum  chemical  study  of  the  binding  of  Li+  cation  to  polyalkyloxides  has  been  carried  out.  The  lithium  cation  interaction  with  three 
poly  alky  loxides  (polyethylene  oxide  (PEO),  polytrimethylene  oxide  (PTMO),  and  polypropylene  oxide  (PPO))  has  been  investigated  using  ab 
initio  molecular  orbital  theory  at  the  HF/6-31G*  level  with  molecular  models  for  the  polymers.  Coordination  by  one  to  six  oxygens  was 
considered.  In  addition,  higher  level  calculations  were  carried  out  using  G3(MP2)  theory  for  coordination  of  Li+  by  one  oxygen.  For 
coordination  of  lithium  by  one  oxygen,  the  binding  energy  ordering  is  PTMO  >  PPO  >  PEO,  with  PTMO  having  the  largest  lithium  cation 
affinity.  The  same  ordering  is  found  for  larger  coordination  numbers  with  the  exception  of  coordination  by  six  oxygens,  where  the  ordering 
changes  due  to  the  steric  interactions. 
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1.  Introduction 

There  has  been  much  interest  in  lithium  polymer  electro¬ 
lytes  for  their  potential  applications  in  fuel  cells,  secondary 
battery  systems,  and  other  electrochemical  devices.  Polymer 
electrolytes  [1,2]  are  generally  composites  of  a  polyethylene 
oxide  (PEO)  or  another  modified  polyether  and  a  salt  such  as 
LiCF3S03  or  Li(CF3S02)2N.  The  ion-polymer  and  ion-ion 
interactions  in  these  materials  play  an  important  role  in  the 
ion  transport  mechanism.  There  have  been  theoretical  stu¬ 
dies  of  lithium  cation-polymer  interactions  [3-15]  and  most 
of  them  have  focused  on  PEO.  We  have  reported  studies 
based  on  a  molecular  model  of  CH3(OCH2CH2)nOCH3  for 
PEO  [3,4  .  The  results  indicated  that  the  binding  energies  of 
Li+-PEO  complexes  increase  with  increasing  coordination 
of  Li+  by  oxygen,  but  the  average  binding  energy  per  Li-0 
bond  decreases.  In  this  paper,  we  report  on  a  comparative 
study  of  the  interaction  of  lithium  cation  with  two  other 
polyalkyloxides  for  comparison  with  polyethylene  oxide. 
The  two  polyalkyloxides  are  polytrimethylene  oxide 
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(PTMO)  and  polypropylene  oxide  (PPO).  Coordination  by 
one  to  six  oxygens  was  considered.  The  lithium  cation 
binding  energies  of  these  polymers  are  of  interest  in  finding 
polymer  electrolytes  with  higher  ionic  conductivities  and  to 
help  understand  the  conductivity  mechanism.  It  has  been 
speculated  that  smaller  binding  energies  may  result  in 
smaller  lithium  cation  migration  barriers,  and  thus,  higher 
ionic  conductivities  [16].  In  Section  2,  the  theoretical  meth¬ 
ods  are  described.  In  Section  3,  the  lithium  cation  binding 
energies  for  a  series  of  polyalkyloxide-Li+  complexes  are 
presented,  and  the  results  are  discussed. 

2.  Theoretical  methods 

All  of  the  chemical  species  considered  in  this  paper  were 
fully  optimized  at  the  HF/6-31G*  level  of  theory  [17]. 
Vibrational  frequencies  were  calculated  to  verify  that  the 
structures  corresponded  to  local  minima.  The  lithium  cation 
binding  energies  of  the  complexes  are  defined  relative  to  the 
corresponding  all-trans  neutral  molecule.  However,  in  some 
cases  other  conformers  have  lower  energies;  and  this  is  noted 
in  those  cases.  The  single  oxygen  coordination  complex  was 
further  studied  at  higher  levels  of  theory  including  use  of 
larger  basis  sets  (6-3 1+G*,  6-31 1G  (2df,  p),  G3MP2Large 
[8])  at  the  second-order  perturbation  level  (MP2)  and 
G2(MP2)  theory  [8]  to  assess  the  effect  on  the  binding 
energies.  All  of  the  calculations  were  done  with  the  Gaussian 
98  [19]  computer  program. 
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3.  Results  and  discussion 

3.1.  Polymer  models 

PTMO  is  modeled  by  CH3(OCH2CH2CH2)nOCH3 
in  =  2-5),  PEO  is  modeled  by  CH3(OCH2CH2)„OCH3 
(n  =  2-5),  and  PPO  is  modeled  by  CH3(OCH(CH3)CH2)^- 
OCH3  ( n  =  2-5).  The  structures  of  the  n  —  2  models  are 
shown  in  Fig.  1.  The  n  =  2  species  are  used  for  lithium 
cation  coordination  with  one,  two,  and  three  oxygens,  the 
n  —  3  species  for  four-coordinated  Li+,  the  n  =  4  model 
for  five-coordinated  Li+,  and  the  n  =  5  model  for  six- 
coordinated  Li+.  The  diglyme  model  ( n  =  2)  for  PEO 
was  previously  [3,4]  used  to  model  one-,  two-,  and  three- 
coordinated  Li. 

Viscometric  studies  [20-22]  show  that  PTMO  is  unu¬ 
sually  compact  presumably  due  to  a  preference  for  the 
gauche  conformation  about  CC  bonds,  leading  to  small 
random  coil  dimensions  (conformational  randomness). 
From  our  calculations,  the  n  =  2  model  for  PTMO  rotation 
about  a  CC  bond  to  form  a  gauche  conformation  results  in 
an  increased  favorable  interaction  between  an  oxygen  atom 
(partial  negative  charge)  and  a  CH2  group  (partial  positive 
charge)  which  are  separated  by  three  bonds.  On  the  other 
hand,  rotation  about  a  CO  bond  results  in  the  unfavorable 


interaction  between  two  large  partially  positively  charged 
CH2  groups.  This  ttgg  conformation,  trans  about  CO  bonds 
and  gauche  about  CC  bonds,  is  responsible  for  the  con¬ 
formational  randomness  of  PTMO  when  averaged  over 
the  repeat  unit.  At  the  HF/6-31G*  level,  the  ttggttgg  con¬ 
formation  of  CH3(OCH2CH2CH2)2OCH3  is  3.2  kcal/mol 
lower  in  energy  than  the  all-trans  species.  The  optimized 
ttggttgg  structure  contains  four  OCH2  interactions  so  each 
gauche  interaction  contributes  about  0.8  kcal/mol  to  the 
enhanced  stability  of  the  structure.  In  the  following  results, 
the  binding  energies  for  the  Li+-PTMO  complexes  were 
calculated  by  comparing  the  energy  of  the  complex  with 
the  energies  of  the  corresponding  all-trans  model  for 
PTMO  and  Li+.  The  all-trans  conformer  was  used  in  the 
binding  energy  calculations  because  it  would  be  extremely 
difficult  to  locate  the  most  stable  conformer  for  the  dif¬ 
ferent  PTMO  models.  Similarly,  the  binding  energies  for 
the  Li+-PEO  and  Li+-PPO  complexes  were  calculated  by 
comparing  the  energy  of  the  complex  with  the  energies  of 
the  corresponding  all-trans  model  for  PEO  and  PPO.  As  in 
the  case  of  PTMO,  there  are  gauche  containing  conformers 
of  the  n  =  2  molecular  models  of  PEO  and  PPO  that  are 
lower  in  energy  than  the  all-trans  conformer,  although  the 
energy  difference  is  much  smaller  (less  than  1  kcal/mol) 
[23-25]. 


PEO 


PTMO 

Fig.  1.  Structures  of  CH3(0CH2CH2)20CH3,  CH3(OCH(CH3)CH2)2OCH3  and  CH3(OCH2CH2CH2)2OCH3  used  as  models  for  PEO,  PPO,  and  PTMO, 
respectively.  The  black  atoms  are  carbon,  the  white  atoms  are  oxygen,  and  the  small  shaded  atoms  are  hydrogen. 
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Fig.  2.  Structures  of  one-,  two-,  three-,  four-,  five-,  and  six-coordinated  Li+-PTMO  complexes.  The  black  atoms  are  carbon,  the  white  atoms  are  oxygen,  the 
large  shaded  atoms  are  lithium,  and  the  small  shaded  atoms  are  hydrogen. 


3.2.  Binding  energies  of  the  Li^-PTMO,  Li+-PEO, 
and  Li^ -PPO  complexes 

One-,  two-,  three-,  four-,  five-,  and  six-coordinated  com¬ 
plexes  of  Li+  with  PTMO  and  PEO  were  investigated.  While 
we  have  not  done  an  exhaustive  search  of  the  potential 
energy  surface  for  conformers,  we  have  tried  to  choose  the 
conformers  that  are  likely  to  be  lowest  in  energy.  The 
geometries  of  these  structures  were  fully  optimized.  The 
structures  of  the  PTMO-Li+  complexes  are  shown  in  Fig.  2. 
The  structures  of  the  PEO  and  PPO  complexes  are  similar. 

The  HF/6-31G*  binding  energies  of  the  PTMO-Li+, 
PEO-Li+,  and  PPO-Li+  complexes  ( n  =  2-5)  are  given 
in  Table  1 .  The  binding  energies  given  in  Table  1  are  relative 
to  the  all-trans  conformer.  The  Li-0  distances  are  given  in 
Table  2.  For  coordination  of  the  lithium  cation  with  a  single 
oxygen  in  the  polymer  model  the  ordering  of  the  binding 


Table  1 

HF/6-31G*  binding  energies  (kcal/mol)  Li+-PTMO,  Li+-PTMO,  and 
Li+-PTMO  complexes21 


n 

Oxygen 

coordination 

Structure 

PEO 

PPO 

PTMO 

2 

1 

1 

39.4 

42.9 

48.6 

2 

2 

66.0 

68.3 

75.0 

3 

3 

87.1 

89.4 

97.4 

3 

4 

4 

103.0 

103.0 

110.4 

4 

5 

5 

110.1 

110.2 

115.6 

5 

6 

6 

115.4 

115.3 

112.7 

a  Using  CH3(OCH2CH2CH2)„OCH3  (n  =  2-5)  as  a  model  for  PTMO, 
CH3(OCH2CH2)/3OCH3(n  =  2-5)  as  a  model  for  PEO,  CH3(OCHCH3- 
CH2)2OCH3  0 n  =  2-5)  as  a  model  for  PPO.  The  binding  energy  is  the 
energy  of  the  complex  relative  to  the  Li+  cation  and  the  all-trans 
conformer.  It  does  not  include  zero-point  energies. 
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Table  2 

Summary  of  Li-O  bond  distances  (in  Angstroms)  in  PEO-Li+,  PPO-Li+, 
and  PTMO-Li+complexes 


n 

Coordination  PEO 

number 

PPO 

PTMO 

2 

1 

1.835 

1.831 

1.830 

2 

1.876 

1.868 

1.851 

1.880 

1.887 

1.860 

3 

1.918 

1.912 

1.876 

1.938 

1.919 

1.897 

1.938 

1.923 

1.897 

3 

4 

2.006 

1.952 

1.937 

2.007 

2.001 

1.954 

2.014 

2.036 

1.968 

2.014 

2.049 

2.004 

4 

5 

2.069 

2.062 

2.042 

2.071 

2.076 

2.047 

2.110 

2.086 

2.093 

2.130 

2.109 

2.103 

2.175 

2.109 

2.181 

5 

6 

2.089 

2.091 

2.074 

2.154 

2.127 

2.074 

2.162 

2.135 

2.207 

2.215 

2.287 

2.208 

2.227 

2.316 

2.426 

2.233 

2.317 

2.426 

energies  is  PTMO  >  PPO  >  PEO.  The  Li+-PTMO  binding 
energy  of  48.6  kcal/mol  is  9.2  kcal/mol  larger  than  that  of 
Li+-PEO.  The  Li+-PPO  binding  energy  is  42.9  kcal/mol. 
The  ordering  does  not  change  if  the  ttggttgg  conformer  is 
used  for  the  n  =  2  model  of  PTMO.  This  conformer  is 
3.2  kcal/mol  lower  in  energy  than  the  all-trans  conformer. 

The  binding  energies  of  the  singly-coordinated  oxygen 
complex  at  several  higher  levels  of  theory  (MP2/6-3 1+G(d), 
MP2(6-31  lG(2df,  p),  G3(MP2)  theory)  are  given  in  Table  3. 
The  results  at  the  higher  levels  give  binding  energies  similar 
in  magnitude  and  order  as  HF/6-31G*.  The  binding  energies 
differ  by  only  several  kcal/mol.  The  stronger  Li+  binding  in 
PTMO  compared  to  PEO  can  be  partially  be  rationalized  by 
examination  of  the  electron  distributions.  The  atomic 
charges  (Mulliken  populations  with  hydrogens  summed  into 


Table  3 

Binding  energies  (in  kcal/mol)  of  PEO-Li+,  PPO-Li+,  and  PTMO-Li+ 
complexes  at  different  levels  of  theory 


Method 

PEO 

PPO 

PTMO 

HF/6-31G(d)a 

39.4 

42.9 

48.6 

MP2/6-31+G(d)a 

38.5 

42.9 

47.2 

MP2/6-311G(2df,  p)a 

41.9 

45.8 

49.8 

MP2/G3MP2Largeab 

37.7 

42.3 

46.5 

G3(MP2)  w/o  ZPEC 

37.9 

42.7 

46.6 

G3(MP2)w  ZPEd 

37.1 

41.6 

45.4 

a  At  the  HF/6-31G(d)  geometries,  no  zero-point  energies  included. 
b  This  basis  set  is  equivalent  to  6-31 1++G(2df,  2p)  for  first-row  atoms. 
See  [18]. 

c  G3(MP2)  theory  without  zero-point  energies. 
d  G3(MP2)  theory  with  zero-point  energies. 


PEO-Li+  complexes  as  a  function  of  coordination  number. 

heavy  atoms)  for  the  single-coordination  structures  shows 
that  the  CH2  groups  two  bonds  away  from  the  central  oxygen 
have  very  different  values  in  the  PTMO  and  PEO  complexes 
(PTMO  =  —  0.04p,  PEO  =  0.30p).  Assuming  a  system  of 
point  charges  with  hydrogens  summed  into  heavy  atoms  (at 
the  optimized  geometries),  the  Coulombic  interaction  ener¬ 
gies  (< qq/r )  between  Li+  and  the  model  polymer  are 
—33.2  kcal/mol  for  the  one-coordinated  PTMO  complex 
versus  —15.5  kcal/mol  for  the  corresponding  PEO  complex. 
Although  the  Coulombic  interaction  energies  are  different 
from  the  HF/6-31G*  results  (—48.6  and  —39.4  kcal/mol, 
respectively),  the  former  are  consistent  with  the  increased 
binding  of  PTMO-Li+.  This  is  primarily  due  to  the  charge 
redistribution  from  the  extra  CH2  group  in  the  PTMO  repeat 
unit.  The  atomic  charges  of  the  PPO  complex  are  very 
similar  to  that  of  PEO  and,  thus,  their  binding  energies 
are  similar. 

A  plot  of  the  binding  energies  in  Table  1  versus  coordina¬ 
tion  number  is  given  in  Fig.  3.  The  plot  indicates  that  the 
binding  energies  of  the  PTMO-Li+  complexes  level  off 
between  four-  and  five-coordination  whereas  they  level 
off  at  around  five-  and  six-coordination  for  PEO-Li+  and 
PEO-Li+.  The  only  six-coordination  structure  that  we  could 
locate  for  PTMO-Li+  had  a  smaller  binding  energy  than  the 
five-coordination  complex.  The  large  range  of  Li-0  dis¬ 
tances  (2.07-2.43  A)  for  six-coordination  PTMO-Li+  sug¬ 
gests  that  the  PTMO  chain  is  having  some  difficulty  in 
forming  a  six-coordination  complex  due  to  steric  crowding. 
The  binding  energies  of  the  PPO-Li+  complexes  are  very 
similar  to  that  of  the  PEO-Li+  complexes.  This  is  expected 
because  of  the  similar  atomic  charges  on  the  models  repre¬ 
senting  PEO  and  PPO. 

3.3.  Conclusions 

In  this  paper,  we  have  reported  on  an  ab  initio  molecular 
orbital  study  of  the  interaction  of  a  single  Li+  cation  with 
molecules  modeling  polyalkyloxides  including  PEO, 
PTMO,  and  PPO.  Coordination  by  one  to  six  oxygens 
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was  considered.  For  the  singly-oxygen  coordinated  structure 
at  the  HF/6-31G(d)  level,  the  PTMO-Li+  complex  had  the 
largest  binding  energy  followed  by  PPO-Li+  and  then  PEO- 
Li+.  This  order  was  confirmed  by  high  level  calculations 
including  G3(MP2)  theory.  Similarly,  for  two-  through  five- 
coordination,  the  binding  energies  of  the  PTMO-Li+  com¬ 
plexes  are  larger  than  those  of  the  corresponding  PEO 
complexes.  This  is  due  to  the  differences  in  the  charge 
distributions  in  the  two  systems.  Only  the  six-coordination 
PEO-Li+  complex  had  a  larger  binding  energy  than  the 
corresponding  PTMO-Li+  complex  due  to  steric  hindrance 
in  the  latter  case.  The  binding  energies  of  the  PPO-Li+ 
complexes  for  the  two-  through  six-coordination  structures 
are  very  similar  to  those  for  the  PEO-Li+  complexes.  A 
comparison  of  Li+  migration  barriers  for  these  complexes 
will  be  reported  on  in  a  future  publication.  A  previous  study 
of  Li+  migration  in  Li+-PEO  complexes  found  that  the 
barriers  largely  depended  on  the  energy  differences  between 
coordination  structures  [4].  Based  on  the  energy  differences 
in  Table  1  the  migration  barriers  for  PTMO  are  likely  to  be 
smaller  than  that  of  PEO  and  those  of  PPO  are  likely  to  be 
very  similar  to  those  of  PEO. 
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